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Abstract 

We consider background dynamics of generalized Galileon theories in the context of 
inflation, where gravity and inflaton are non-minimally coupled to each other. In the 
inflaton oscillation regime, the Hubble parameter and energy density oscillate violently 
in many cases, in contrast to the Einstein gravity with minimally coupled inflaton. 
However, we find that there is an adiabatic invariant in the inflaton oscillation regime 
in any generalized Galileon theory. This adiabatic invariant is useful in estimating 
the expansion law of the universe and also the particle production rate due to the 
oscillation of the Hubble parameter. 
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1 Introduction 


Inflation [T] is now widely accepted as a successful way of producing the flat and isotropic 
universe as well as small density perturbations as seeds for galaxies. It is triggered by the 
potential energy of a scalar held called inhaton, which slowly rolls down its potential during 
inhation [2]. For successful reheating, inhaton inevitably couples to light particles, and this 
coupling may induce sizable radiative corrections to its ehective potential. Even if this is not 
the case, there might be a self interaction term of the inhaton. As a result, it might contradict 
with the present observational data by producing too large density perturbations [3|. Among 
several approaches to this issue, one way is to generalize the interaction between inhaton 
and gravity to hatten its ehective potential. However, the coupling of inhaton to gravity 
in general produces higher derivatives in their equations of motion. Since it can induce 
other degrees of freedom which might lead to ghost instabilities, it would be desirable to 
construct models without introducing new degrees of freedom. There is a class of theories 
with a scalar held and gravity, called generalized Galileons |1], which do not suher from such 
higher derivatives. Ehorts have been made to study general setups for using these theories 
in inhation (generalized G-inhation [5]), and such theories are under close investigation 
especially after the discovery of the Higgs boson at the LHG PEI and after the results from 
Planck satellite P, in the context of identifying the inhaton as the Higgs held [SHT^. 

In this paper, we study the oscillation regime of the generalized Galileon theories in the 
context of inhation, because it is necessarily in order to understand the full reheating dynam¬ 
ics and construct a consistent cosmological scenario. In fact, it has recently been pointed out 
that in some models of the generalized G-inhation the Hubble parameter oscillates violently 
during the inhaton oscillation regime [13] , which makes the analysis of expansion law of the 
universe non-trivial and also it may lead to Laplacian instability [H] . The knowledge of the 
expansion law of the universe in these models is also necessarily to make a precise prediction 
for the spectral index of the density perturbation, for example. For these purposes, we con¬ 
struct a useful quantity, which we call an adiabatic invariant, in analyzing the background 
dynamics. The expansion law of the universe during inhaton oscillation regime is easily 
derived by using this invariant. It is also useful to extract the oscillation part of the Hubble 
parameter, which helps us estimate the particle production rate via gravitational ehect. 

The organization of this paper is as follows. In Sec. El we explain the procedure to 
obtain the adiabatic invariant, to express the oscillation of the Hubble parameter in terms 
of the inhaton, and to obtain the expansion law of the universe. In Sec. [31 we illustrate this 
procedure using some examples. The last section is devoted to conclusion. 


2 General discussion 

We consider the following action throughout this paper: 

S = Sg + Sm, ( 2 . 1 ) 
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where Sg is the action for the gravity and the inflaton (j), including non-minimal couplings 
between them, and Sm is the action for the matter heldl^ We assume that Sg is given by 
the Lagrangian for generalized Galileon theories [1]: 

r _ ® 

Sg= (2.2) 

d i=2 

where 

C2 = G2{(j),X), (2.3) 

/:3 = -G'3(0,x)n0, (2.4) 

A = ^4(0, X)R + G,x [(□</>)' - (V^V,0)2] , (2.5) 

£5 = G5(0,X)G'^"V^V.0- ^ [(□0)=' - 3(n0)(V^V.0)2 + 2(V^V.0)3] . (2.6) 

Here Gj’s are general functions of 0 and X = —5f^^V^0Vi/0/2, and the subscript X de¬ 
notes the derivative respect to X. Moreover, (V^V,y0)^ and (V^Vi,0)^ are understood as 
V^Vi.0V^V^0 and V^V,y0V'^V^0VpV^0, respectively. We take the background part of 
the metric to be the Friedmann-Lemaitre-Robertson-Walker (FLRW) one with negligible 
curvature, 

= —dt^ + a{tY5ijdx'^dx\ (2.7) 

where a is the scale factor. We also assume that the background part of the inflaton is 

homogeneous, 7 e. 0 = 0(t). In such a case, the background parts of the actions for Eqs. fl2.3p - 
fl2.6l] are respectively given by 


£2 — G 2 , (2-8) 

£3 = —£^300^ ~ £^3X0^0, (2.9) 

£4 = - GH^Gi - 6i7£'400 + 6i72£'4x0^ (2.10) 

£5 = -3772G5<a02 + (2.11) 


after integration by partJ^. where H = d/a is the Hubble parameter, the dot denotes the 
derivative with respect to t and the subscripts of X and 0 are the derivative with respect 
to those quantities. Note that all terms except for £3 depend only on H, 0 and 0. In fact, 
the 0 dependence of £3 can also be eliminated by integration by parts0 Thus, we take the 


Though we focus on the oscillation regime of single-field inflation, generalization of the adiabatic in¬ 
variant to multi-field case is trivial since, as we will see later, the necessary equation is only Eq. (12.161) . 

It is known that some of these models have the ghost instability after inflation, and thus it is even 
non-trivial whether the inflaton can oscillate around the minimum of the potential after inflation |15] . The 
terms which potentially lead to such ghost instabilities are understood as those with 0^ in Eqs. (|2.8I) - (I2.11I) 
whose coefficients are not positive definite. In this paper, we simply neglect such models because they may 
not be connected to the standard big-bang cosmology. 

This can be done for G 3 polynomial for 0 and X as 


0'"X”002 


1 

2n + 3 



1 

2n -I- 3 



3i70’”02”+3 -k m0’”-i0^”+'‘ 


( 2 . 12 ) 
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background part of So generally as follows: 


Sg = 


d'^x a^C{H, <j), (j)). 


(2.13) 


In this paper, we simply neglect the back reaction to the inflaton due to particle produc¬ 
tion]^ We also assume that £ is a sum of polynomials with respect to its arguments. 

The background equations of motion for the system fl2.13l) is given by 


(£^)' + 3HC^ - = 0, 

(2.14) 

o' 

J 

1 

1 

(2.16) 

{Ch)' + 3HCh -3C = 0, 

(2.16) 


where the subscript denotes the derivative with respect to that quantity. The first equation 
is 0’s equation of motion and the second one is the Friedmann equation. In addition, the last 
equation is obtained by taking functional derivative with respect to the scale factor. Note 
that, as we will see below, the last equation is useful in analyzing the inflaton oscillation 
regime although it is redundant. 

In this paper, we focus on the inflaton oscillation regime. We define rries as a typical 
frequency of the inflaton oscillation in the following discussion. In general, nies may be 
estimated as mes ~ 10/01- In the inflaton oscillation regime, the inequality mes > H is 
satisfied. As we will explicitly see in the next section, the Hubble parameter and the energy 
density of the inflaton violently oscillate with time in some models, i.e. H = 0{mesH)- 
However, there is a quantity which evolves only with a time scale of H~^, not m~g, for the 
general Lagrangian fl2.13|) . Such a quantity, which we call an “adiabatic invariant,” can be 
used to evaluate the cosmic expansion law in the oscillation regime. In addition, it is also 
helpful in estimating the amount of particle production via gravitational effects. 


2.1 Adiabatic invariant 

We define an adiabatic invariant Q by the following equality: 


Q = 0{HQ), 


(2.17) 


which implies that Q evolves only with the time scale of From the Friedmann equation 


fl2.15p . we can see that C ~ HCh is satisfiedl^. Then, from Eq. 02.161) we obtain 


{Ch)' = 0{HCh). (2.18) 

It is known that some of these models have a Laplacian instability coming from the negative sound 
speed squared for the scalar perturbation [T4l[T5], and such an instability may lead to a rapid decay of the 
inflaton condensation. Anyway, in order to determine whether there is a Laplacian instability or not, we 
must first solve the background dynamics. 

In fact, if £ ^ HCh does not hold, C and must cancel almost exactly with each other. This occurs 
only when £ is linear in (j), and we do not consider such a case here. 
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Thus we have proven that Ch is an adiabatic invariant. This invariant can also be expressed 
by 0 and 0, by eliminating H using Eq. 02.151) 

Here we dehne 


#6 


J = 


c 


H 


ml' 


(2.19) 


and 


pj = 3 MIj‘^ 


( 2 . 20 ) 


Both J and pj are adiabatic invariants. Note that they coincide respectively with the Hubble 
parameter and the energy density in the case of the minimal Einstein gravity. In this sense, J 
and Pj may be interpreted as generalization of the Hubble parameter and the energy density 
to the case of the generalized Galileon theories, respectively. 


2.2 Evaluation of the cosmic expansion law 

In this subsection, we schematically describe the procedure to evaluate the cosmic expansion 
law by using the adiabatic invariant that we have derived just above. Some specihc examples 
of this procedure will be given in the next section. 

In order to obtain the expansion law, we hrst evaluate the scale factor dependence of 
Ch- Let us start with Eq. fl2.16j) . We have to express the oscillation average of C by that of 
HCh- The oscillation average of Eq. 02.151) gives 

{£) - - {HCh) = 0, (2.21) 

where the angular bracket denotes an oscillation average which smears out the oscillation 
with the frequency of the order of ~ rries- In addition, by multiplying Eq. 02.14p with 0 and 
taking an oscillation average, we obtain the Virial theorem: 

(.#.£0 + ~ 0 , ( 2 . 22 ) 

where we have taken only leading terms with respect to 0{H/meg). Here, note that usually 
C can be decomposed into at most three terms which are proportional to (pC^, (pC^ and 
HCh, if one neglects small terms in each regime under consideration. See also examples 
in Sec. [3l Thus, from Eqs. 02.211) and 02.22p together with such a decomposition, one can 
express (£) in terms of {HCh)- Then, Eq. 02.16p can be rewritten by taking the oscillation 
average as 

j + c(Tr)J = 0, (2.23) 

where c is 0{1) numerical coefficient which can be determined explicitly if we specify the 
model. Here we used the fact that J is an adiabatic invariant, and thus {HJ) = {H) J. Also 
remember that J is proportional to Ch- Eq. 02.23p can be used to express J as a function of 
the scale factor a. Since we know the relations among {HCh), {4’^<i>) and {(f>C^), in general 
we can express J only by (if), i.e. J (x {H)"^, where m is some number. Then, the scale 
factor dependence of {H) and hence the expansion law of the universe can be estimated. We 
will see how this procedure works for some specihc models in the next section. 

Note that the scale factor a is also an adiabatic invariant by definition. 
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2.3 Particle production 

As we said repeatedly, the Hubble parameter can violently oscillate with time for the general 
Lagrangian fl2.13p . In such a case, the violent oscillation of the background metric causes 
production of non-Weyl invariant particles [16]. Here, we briefly summarize basic ingredients 
of such a particle production event. More detailed discussion is given in na 
We call the produced (scalar) particle y here. The parameter, q, is dehned as 


q = 



(2.24) 


where we have assumed that the v’s mass squared has an oscillating part Am^ due to the 

X, 

inflaton oscillation. Then, the decay rate of 0 is estimated as 




$2 

c 


(2.25) 


where <hc is the canonically normalized amplitude of 0. In the case where y is minimally 
coupled to gravity and the bare mass of y is negligible, the effective mass of canonically 
normalized y is given by the Hubble parameter squared and its time derivative since the 
normalization of y depends on the scale factor. Hence, the decay rate is given by 

{sny 

~ mail' 

where 6H is the oscillating part of the Hubble parameter. Here, we have used the fact that 
6H > is always satished in the inflaton oscillation regime. 

Note that by using the fact that J is an adiabatic invariant, the oscillating part of the 
Hubble parameter can easily be deduced from the relation between J and H. This is another 
beneht of the adiabatic invariant. We will see this point for some specihc models in the next 
section. 


(2.26) 


3 Examples 

In this section, we show several examples to illustrate the procedure explained in Sec. [2l 
We assume that the inflaton potential (the X-independent part of —G 2 ) has the following 
monomial formic 

V(<f) = -r, (3.1) 

n 

and that the inflaton is oscillating around the minimum of the potential with amplitude <h. 

Even for a more general form of potential, we can approximate the potential as the dominant part of 
it. Then, the situation typically reduces to the monomial potential 
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3.1 Example 1 

We first consider 


G2 = W-i/(0), 

G4 = /( 0 ), 


where /(0) is assumed to be of the form 


/(0) = M| 1 + A 


Mp 


(3.2) 

(3.3) 


(3.4) 


Here 0 is understood as the deviation from the potential minimuni^^l. In this model 0 starts 
oscillating around 0 ~ Mp//i, therefore we assume |/i0|/Mp <C 1. 


3.1.1 Adiabatic invariant 

The background action is given by 


Sg= a? 


( 1 + /i 


Mf 




after integration by parts. Thus, Ch is given by 


Ch = 1 + A 


Mf 


- ^Mlh 


Mp' 


and hence the adiabatic invariant J is 

j = w (1 + A 


0 \ , /l 0 


Mf 


+ 


2 Mf 


(3.5) 


(3.6) 


(3.7) 


3.1.2 Expansion law 

Since we assumed |A0|/^p "C 1, the Lagrangian is approximated as 


C - -3MlH^ - 3MlHh^ + ^02 _ (3 8 ) 

IVIp z 

Thus, it can be decomposed as 

Strictly speaking, 0 = 0 does not correspond to the minimum energy for the homogeneous 0- 
condensation, around which 0 oscillates, since the stationary point of 0 receives contribution from fi(j)/Mp 
term. However we neglect such a small effect here. 
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Figure 1: Time evolution of H and J for Example 1 is shown. In these figures, we have taken 
V = Ai?!)^/4, A = 1 and /i = 0.2. H oscillates violently {Left) while J does not {Right). 


Note that this equation is the same as the one in the minimal case of Einstein gravity, 
where the second term in Eq. 03.81) does not exist. Eqs. 02.211) and 02.22p together with the 
decomposition Eq. 03. 9 p gives 




(3.10) 

Substituting this into oscillation-averaeed Ea. 02.1bD srives 




(3.11) 

and hence 

we obtain 



3n 

J (X a "+2 . 

(3.12) 

Note that 

= 0 , and thus the expansion law of the universe is estimated as 



3n 

{H) J X a "+2 , 

(3.13) 

or 

(//).”+ T 

' ^ 3n t 

(3.14) 


This expansion law is the same as that in the minimal case of Einstein gravity. We plot the 
time dependence of H and J (multiplied by f) for n = 4 in Fig. [TJ 


3.1.3 Particle production 

Recalling that J is conserved during the inflaton oscillation, the oscillating part of the Hubble 
parameter is deduced from Eq. 03.7p as 


6H ~ 



(3.15) 































Then, from Eq. fl2.26jl . the production rate is estimated as 


r 




ft 


2 ^eff 


where rueff ~ This reproduces the result in [min]. 

3.2 Example 2 

As a second example, we consider 

G2=X-V{(j)), 

Ga = ImI+ ^ 


^ 2 '''^ ' 2M2' 


Note that this system is equivalent to 


G,=X- V(,p), 

G, = hi}., 

= 


2 M 2 ’ 

where Mp is the reduced Planck mass and M is some mass parameter. 

3.2.1 Adiabatic invariant 

The background action becomes 


Sg= (f^x a? 


+ ( 1 + 


2\ J,2 


3H^ 

Jp I 2 


- 1 / 


after integration by parts. Therefore, Cp is given by 


Ch = -QMiH 1 - 


2M|,M2 / ’ 


and the adiabatic invariant J is given by 

J = H { 1 - 


2M2M2 ] 


(3.16) 


(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 


(3.22) 


(3.23) 


(3.24) 
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3.2.2 Expansion law 

The background dynamics of this system is extensively studied in [T^IT^ . where it is shown 
that this system undergoes a nontrivial regime when S> 1. Let us focus on such a 

regime, when the Lagrangian is approximated as 

C ~ - V. (3.25) 

Then, it can be decomposed as 

£ ~ (3.26) 

Eqs. fl2.2ip . fl2.22p and fl3.26p give 

Substituting this into oscillation-averaged Eq. fl2.16p . one obtains the evolution equation for 
the adiabatic invariant J: 

• St? 

From this equation one can derive 

3rt 

J (X a 2 (n+i) _ (3.29) 


Now let us estimate the expansion law of the universe. Since J x H{1 — 0^/2MpM^), and 
the amplitude of 0^/MpM^ is constant during the regime we consider we have 


_ on 

{H) cx a 2 (^+ 1 ) ^ 

and hence the expansion law is given by 

2(n+ 1) 1 


{H) = 


3n t 


(3.30) 


(3.31) 


We plot the time dependence of H and J (multiplied by t) for u = 2 in Fig. |2] 


|!:10 


3.2.3 Particle production 

From Eq. 03.241) . the oscillating part of the Hubble parameter is deduced as 


5H ~ 


J02 

MlM^' 


(3.32) 


This can be seen by the Friedmann equation. See [T31II1] for detail. 

The center of the oscillation in the left panel of Fig.[2]is not 1, but time-averaging of H gives {H) t = \. 
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Figure 2: Time evolution of H and J for Example 2 is shown. In these figures, we have taken 
V = nr?(^/2 and m/M = 10^. Left: Plot of Ht. H violently oscillates in the non-minimal regime 
H^IM^ ;§> 1. The expansion law is {H)t = 1 for this regime, while it gradually approaches to 
the value with canonical kinetic term 2/3 for <C 1. Right: Plot of Jt. It has a suppressed 

amplitude of oscillation compared to Ht. 


Thus, the decay rate of 0 is estimated as 






(3.33) 


where we have taken J ~ as an order estimation, and <Fc ~ H^/M. Also, rueff ~ 
This reproduces the result shown in |TTj (for case B, wh ere t he non- 
minimal coupling effect is larger than the one that exists in the minimal setup) 


3.3 Example 3 


Finally let us consider the system 


where M is some mass parameter. 


G2 

G 4 

G 5 


x-y(0), 

(3.34) 

04, 

(3.35) 

0 ^ 

M4’ 

(3.36) 


3.3.1 Adiabatic invariant 


The background action is 


Sg = 


(fx a? 





M4 ) 


2 


1 / 


(3.37) 


In fact, from the Friedmann equation, one can show that /MfM^ ^ t2(l) is always satisfied for 

H ^ M in this model. Then, the decay rate reduces to F^-^^ ^ which is the same as that 

obtained in [Hi¬ 


ll 






















from which Ch is given by 


Ch = -6M^H 


/ 

Mj,M^ j 


Then, the adiabatic invariant J is given by 


J 


H 1 


3(^202 \ 

Mp^ ) 


(3.38) 


(3.39) 


3.3.2 Expansion law 


We focus on the regime when the non-minimal kinetic term takes a dominant role here. 
Then, the Lagrangian is approximated as 


im 


2x2 i2 


c ~ p - V, 


M4 2 


and hence it is decomposed as 


£ ~ 

2 n n ^ 


Therefore, Eqs. fl2.22p . fl2.2ip and fl3.4ip give 


(^) = (HCh) . 


2{n + 2) 


Substituting this into Eq. fl2.16p . we get 


3n 




and thus 


_ 3n 

J (X a 2(n+2) 


(3.40) 

(3.41) 

(3.42) 

(3.43) 

(3.44) 


Since the absolute value of the second term in the parenthesis in Eq. fl3.39p is less than 1/3, 
as is guaranteed by the Friedmann equation, this term does not affect the power-law form 
of J. Thus, 

_ 3rr 

{H) oc a 2 ("+ 2 ) (3.45) 


and the expansion law of the universe is given by 


{H) 


We plot H and J (multiplied by t) for n 


2 (n + 2) 1 
?>n t 


2 in Fig. 


tl2 


(3.46) 


As seen in Fig. [31 there are spikes in the Hubble parameter. These correspond to d — 0, and their 
typical time scale is ^ TOeff/m^ (<C 1/meff). At these spikes the quantity in the parenthesis in the RHS of 
Eq. (13.391) changes from 2/3 to 1 and thus H changes from the maximum to the minimum of its oscillation. 
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Figure 3: Time evolution of H and J for Example 3 is shown. In these figures, we have taken 
V = m/Mp = 1 and M/Mp = Left: Plot of Ht. H violently oscillates in the 

non-minimal regime 3> 1. The expansion law is {H)t = 4/3 for this regime, while it 

gradually approaches to the value with canonical kinetic term 2/3 for <C 1. Right: Plot 

of Jt. It has a suppressed amplitude of oscillation compared to Ht. 


3.3.3 Particle production 

From Eq. fl3.39p . the oscillating part of the Hubble parameter is 




5H ^ J 




(3.47) 


and hence the decay rate is estimated as 


"" Mj,M^ ’ 

where we have taken J ~ 77 as an order estimation, and ~ 

(MV77<h)A^/2^’^/2-i_ 


(3.48) 


Also, rUeff 


4 Conclusion 

In this paper, we study the generalized Galileon theories in the context of inflation, especially 
focusing on the inflaton oscillation regime. During this regime, in general, these models lead 
to violent oscillations of the Hubble parameter and the energy density of inflaton, and hence 
it would be desirable to construct an adiabatic quantity which evolves only with the time 
scale of H~^. We have shown the existence of an adiabatic invariant Q = 0{HQ) during this 
regime, and also shown an explicit method to construct this adiabatic invariant. To conhrm 
our procedure, we have demonstrated that the proposed adiabatic invariant actually evolves 
with the time scale of in several examples of generalized Galileon theories by numerical 
calculation. Such an invariant is useful in estimating the expansion law of the universe, and 
expressing the oscillation in the Hubble parameter in terms of the inflaton so as to estimate 

#13 production rate may be enhanced by the effect of spikes mentioned in the previous footnote. 
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the amount of particle production via gravitational effects. Thus it may take crucial roles 
to predict the precise value of the spectral index. 

We also emphasize that the analysis background dynamics is a hrst step to understand 
the actual phenomena in the oscillation regime. In a specihc example of G^'^d^(j)dy(j) model, 
the relation \H\ ~ mesH ^ was crucial for the Laplacian instability ng. Without 
correct understandings of the background dynamics, we cannot judge the phenomenological 
validity of various inflation models. 
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